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Abstract 

We introduce an improved version of the simulation code FEWZ (Fully Exclusive W 
and Z Production) for hadron collider production of lepton pairs through the Drell- 
Yan process at next-to- next-to-leading-order (NNLO) in the strong coupling constant. 
The program is fully differential in the phase space of leptons and additional hadronic 
radiation. The new version offers users significantly more options for customization. 
FEWZ now bins multiple, user-selectable histograms during a single run, and produces 
parton distribution function (PDF) errors automatically. It also features a signifcantly 
improved integration routine, and can take advantage of multiple processor cores locally 
or on the Condor distributed computing system. We illustrate the new features of 
FEWZ by presenting numerous phenomenological results for LHC physics. We compare 
NNLO QCD with initial ATLAS and CMS results, and discuss in detail the effects of 
detector acceptance on the measurement of angular quantities associated with Z-boson 
production. We address the issue of technical precision in the presence of severe phase- 
space cuts. 



1 Introduction 



Electroweak gauge boson production is a standard candle for Large Hadron Collider (LHC) 
physics studies. It is a background to Z' production and numerous other new physics 
searches. As one of the cleanest processes with copious production (millions per year at 
design luminosity), it can be used as a luminosity monitor [lj, to constrain PDFs (parton 
distribution functions) [2], and to study electroweak physics parameters [3]. Therefore, un- 
derstanding its production is crucial as the LHC physics program moves forward. With so 
many events, systematic errors dominate statistical ones, and are expected to eventually 
reach 1 — 2% [3] at the LHC. Next-to-leading order (NLO) predictions in the strong cou- 
pling, with 0(10%) errors, are insufficient for a precise comparison with data; more accurate 
calculations are required. Predictions through NNLO in perturbative QCD must be used. 

The inclusive 0(a 2 s ) corrections to electroweak gauge boson production have been known 
for some time [5]. The theoretical uncertainties are at the percent level. Exclusive produc- 
tion, which is necessary for any realistic prediction or phenomenological study in a detector 
of finite acceptance, is technically challenging but has been achieved [&HTTj. One of us pre- 
viously released a public simulation code FEWZ (Fully Exclusive W and Z Production) 
that implemented the NNLO predictions and allowed for arbitrary kinematic cuts to be 
imposed. However, this previous version suffered from several shortcomings. It was not 
easily customizable; only one cross section of interest could be calculated at a time, thereby 
necessitating numerous runs of FEWZ to obtain a kinematic histogram. For severe cuts 
on the leptonic phase space, a Monte Carlo integration error below a few percent was not 
achievable [T2~|[T3"]. 

In this manuscript we present a new version of FEWZ which addresses the issues described 
above. Specifically, the new features of FEWZ are listed below. 

• The user can define multiple, arbitrary kinematic variables to be binned automatically 
during a single run. Most of the commonly desired histograms are included in the new 
distribution of FEWZ. 

• The calculation has been broken up into 230 sectors that can be run in parallel, dra- 
matially improving the speed and final numerical integration error. Sub-percent in- 
tegration errors are easily obtainable even in the presence of significant phase-space 
restrictions. 

• For all current PDF sets, errors are automatically calculated for the total cross section 
and all histogram bins. 

• Most parameters of interest, such as cuts and couplings, are now set in an external 
input file, allowing the user complete control over the settings of a run. 

We focus on the production of through 7* ' jZ in this manuscript; W production will be 
addressed elsewhere. To demonstrate the new FEWZ, numerous phenomenological results 
for the LHC are presented. A detailed study of PDF and scale uncertainties is performed 
for numerous observables. We discuss the effects of acceptance cuts and theoretical errors 
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on the measurement of angular quantities in Z-boson events, such as the Collins-Soper 
angles and moments [HHT6] , and a proposed transverse-plane angular cut designed to reduce 
experimental backgrounds. A comparison with initial ATLAS and CMS results is performed. 
We show explicitly that the technical limitations found in Ref. p2] are solved. Further details 
on how to install and run FEWZ are available in the accompanying manual. 

Our manuscript is organized as follows. In Section [21 we review the calculation of the 
NNLO corrections and their implementation in FEWZ. Section [3] describes the new version 
of the code, with emphasis on the improvements. Phenomenological results are presented in 
Section HI We conclude in Section [51 

2 Review of FEWZ 

The version of FEWZ we consider here calculates the fully differential production of dilepton 
pairs via the neutral current (intermediate photons and Z-bosons). It is designed to make 
predictions for hadron-collider observables with realistic acceptance cuts at NNLO in the 
strong coupling constant. All spin correlations and finite- width effects are included. The 
residual scale error on typical cross sections is less than 1% and is in good agreement with 
the NLO scale error band. For more details, we refer the reader to Refs. [HE]. 

2.1 Calculational details 

In QCD factorization, the Drell-Yan differential cross section for hih 2 — >■ 7*,Z — > lil 2 can 
be expressed as 

d& = ^2 dx i dx 2ft 1 {xi)fj 2 {x 2 )dcr ij ^ lll2 (x 1 ,X2), (1) 

ij J 

where the PDFs /j represent the probability of obtaining parton i from hadron h, and da^ 
represents the partonic cross section. The leading order (LO) and NLO contributions are well 
known [T7j. To calculate the NNLO corrections, three types of contributions must be 

considered: two-loop double- virtual contributions, one-loop real-virtual contributions with 
the emission of an extra parton, and tree-level double-real contributions, with emission of 
two extra partons. Each piece is separately divergent, and must be summed to obtain a 
finite result. 

The loop integrals of the first two types, double-virtual and real-virtual, are dealt with 
by decomposing the Feynman integrals into a basis of so-called master integrals in an auto- 
mated fashion [18] . These master integrals can be expressed in dimensional regularization as 
a Laurent series in the dimensional regularization parameter e in terms of known functions. 
The double- real contributions are also divergent, due to soft and collinear singularities. The 
singularities are extracted by mapping each propagator denominator into a set of hypercube 
variables, each of which only vanishes at one endpoint. Since each denominator typically van- 
ishes when multiple variables reach an endpoint, the technique of sector decomposition [19] 
must be used to obtain the required form. This method involves splitting the integrand into 
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multiple terms, called sectors, which correspond to the different singular limits of the pro- 
cess. For more details for the sector decomposition required in this process, see Ref . [20l|2T] . 
The result of this process is a separation of the singular limits which allows them to be 
independently extracted. One can perform the expansion 



where x represents a hypercube variable, for each propagator topology. The end result 
is approximately 200 sectors corresponding to different initial partons, real/virtual pieces, 
soft /collinear counterterms, and mappings of the hypercube variables corresponding to the 
singularity structure of the matrix elements. The coefficients of the pole terms in e are 
numerically checked to cancel when summed, and the code calculates only the e° coefficient 
for each sector. The original calculation is more fully described in jSKH] . We note that FEWZ 
incorporates only QCD corrections. Effects such as final-state photonic radiation which can 
be important for some experimental cuts are not simulated. 

2.2 Details of the numerical code 

The resulting finite contributions are implemented in a Fortran code and interfaced with 
routines for calling PDF sets to complete the hadronic calculation. The list of supported 
PDF sets has been greatly expanded and will be discussed in the next section. This results 
in a (4, 7, 11) dimensional integrand for the (LO, NLO, NNLO) computation, corresponding 
to the allowed degrees of freedom of the observable leptons and jets; there is an extra param- 
eter in the NNLO integrand for internal purposes. This is then interfaced with an adaptive 
Monte Carlo numerical integrator. Because the kinematics of all final-state particles are re- 
constructed in the integrand, cuts can be imposed by checking whether requirements are met 
and zeroing the integrand appropriately. In the previous version, these cuts were input into 
a Fortran source file, requiring recompilation if one wants to run with a different set of cuts. 
Several other parameters of interest, such as masses, couplings, renormalization/factorization 
scale, collision energy, and collision type (proton-proton vs. proton-antiproton), are set in 
an external input file in both the old and new versions of FEWZ. 

For the neutral current, one is often interested in the Z-resonance. In this case, a variable 
transformation 



is made, which has the effect of flattening the Z-resonance, greatly improving the efficiency 
of numerical integration. The new version of FEWZ allows this variable change to be turned 
on or off, depending on the region of lepton invariant mass under consideration. To reach a 
numerical precision better than the errors due to scale variation (1%), the previous version 
of the code typically required days to run for non-trivial cuts on the leptonic phase space. 
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3 Description of code improvements 



We have rewritten and significantly improved FEWZ. The new version features an improved 
numerical integration routine and offers users significantly more options for physics studies. 
The resulting program is a powerful tool that enables precision studies of most aspects of 
lepton-pair production at hadron colliders. We describe in this section the new features of 
FEWZ. More details on running FEWZ can be found in the accompanying manual. 

3.1 Overview of improvements 

Our neutral-current code has been updated to improve speed and precision, as well as increase 
the amount of information generated during a single run. We provide a summary of the 
changes and their intended effect in this section. While we focus here on neutral-current 
production, the same improvements apply also to W production. This will be addressed in 
a future publication. 

Parallelization: Each of the 230 NNLO sectors is calculated independently. This allows 
the Monte Carlo integration to adapt to the structure of a single sector rather than 
to all at once, and also allows more than one processor core to work simultaneously. 
We have written scripts for starting multi-core local runs, and also for running on the 
Condor [22] distributed computing system, as well as combining results from individual 
sectors. 

Run parameters: All inputs, including cuts on leptons and jets, electroweak couplings, 
and other parameters which control run setting, are now set in an external input file, 
allowing the user complete flexibility to customize FEWZ. 

Histograms: By tabulating the weights associated with each event, kinematic distributions 
are now produced automatically during a run, with little overhead. The user can select 
which histograms to fill in an external input file. Most distributions of interest are 
included in the default version of FEWZ. 

PDF errors: When running with PDF sets that contain error eigenvectors, all eigenvec- 
tors are calculated automatically for each histogram bin. The resulting output can 
be combined using the included scripts to produce a final output file that contains 
the integration error as well as PDF error for both the total cross section and each 
histogram bin. 

3.2 Details of the numerical integration 

Quantities of interest, such as the total cross section and its various kinematical distributions, 
are produced by numerically integrating Eq. flTJ with a Monte Carlo adaptive integrator. We 
use the standard Vegas routine from the package Cuba 1.7 [23J, which is distributed with 
our program. Vegas allows one to save the state of the integration between grid adaptations, 
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which is useful for long calculations and allows a pause so that intermediate output is pro- 
duced. In addition, the weight of each sampling point is returned, which we use to calculate 
histograms bins as described later. 

Several additional variable transformations are implemented to improve the performance 
of the numerical integration. In addition to the removal of the Z propagator using Eq. fl3]), 
an additional smoothing of certain sectors' integrands is performed. After the e-singularities 
are removed, the integrands may still diverge logarithmically at or 1 in the hypercube 
variables, even though the integrals are still finite. For stability of the adaptive integrator, 
transformations such as 



are performed for the NNLO sectors. Such a transformation removes singularities of the 
form ln(x) while restricting the integration region's support to the unit hypercube. 

3.3 Parallelization 

As mentioned, the integrand has been broken into 230 sectors corresponding roughly to the 
results of sector decomposition. Some sectors, however, have been merged or split. After 
sector decomposition, there were approximately 260 sectors; these correspond to different 
initial state partons (gluon-gluon, gluon-quark, quark-quark), different diagram type (real- 
real, real- virtual, soft /collinear counterterms) , and different decomposition of the singularity 
structure. A dummy run of these sectors was set up for a typical input. It was found that 
some sectors anti-correlate over the randomly generated phase space points. For the results 
of such sectors such i and j, this implies 



where cr, denotes the standard deviation of the random sample of integrand i, and Cov(i,j) 
the covariance of the two sectors. This indicates that these sectors should be combined, since 
the resulting error would be smaller than random sampling each independently, i.e., they 
are canceling. Candidate combinations were tested with multiple sets of cuts, and those that 
performed better together were combined. Keeping the remainder separate allows Vegas to 
adapt its grid to the different integrand shapes better and accelerates convergence. 

A few sectors typically take much longer to reach a target error than others due to high 
variance and evaluation time. To prevent these from holding up the user in a cluster envi- 
ronment where numerous processors are available, identical copies of these extreme sectors 
are split over multiple new sectors with differing random seeds to statistically reduce their 
contribution to the total error. Each sector is given a target precision to reach depending 
on the total absolute precision set in the input file. In the worst-case where the integrator 
fails to adapt (e» ~ yJVi/N), it is found to be most CPU-efficient to weight the goal error 
for each sector according to 



dx — > 6w(l — u)du 



(4) 




(5) 



where Vi is the sector's random-sampling variance (estimated for default cuts) and is the 
evaluation time per point. 

The sector number can be set in the input file and run with the compiled program, which 
will then produce a human-readable output file. However, it is intended that the user instead 
run a provided starting script which will create a directory structure for the results of all 
sectors and run the code on all sectors. There is a script for running locally on a specific 
number of cores, and one for running on the Condor system, where a job is submitted for 
each sector. The user can then run a combining script once output files have been produced 
for each sector; this only requires that one iteration of Vegas has completed. Further details 
are available in the FEWZ manual. 

3.4 Run parameters 

All relevant input parameters are now set in an external input file. These include vector- 
boson masses and couplings of the photon and Z to fermions. The vector and axial couplings 
of the Z are set separately from cvqed and sin 2 ^, allowing an arbitrary electroweak scheme 
or an improved Born approximation to be implemented. The user may also select whether 
to optimize the integration for the Z-peak as described in Eq. fl3]). The new version of 
FEWZ allows numerous cuts to be selected in the input file. These include restrictions on 
the following quantities: 

• lepton transverse momenta and the dilepton transverse momentum; 

• lepton pseudorapidities and the dilepton rapidity; 

• dilepton invariant mass; 

• jet transverse momenta and pseudorapidities; 

• the number of observable jets; 

• jet-jet, jet-lepton, and lepton-lepton isolation. 

The input file also allows one to select the jet-merging algorithm (cone or anti-Zc-r), as well 
as the chosen PDF set. We detail in Sec. 13.61 the supported PDFs. 

3.5 Histogramming 

Since each point in the 11-dimensional parameter space generated by Vegas corresponds to 
particular kinematics, we can save this information to reconstruct more detailed distributions 
than just the total cross section, with little overhead. The "event" (whose contribution may 
be negative) is sorted into the appropriate bin for each of the many distributions defined. 
The bin size and histogram extent can be changed by the user in a histogram input file. 
Since Vegas returns the weight of the generated point, this weight is used to keep track of 
the weighted average and standard deviation for each bin in the same fashion as Vegas for 
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the total. A copy of each histogram bin is made for each PDF error eigenvector, reweighted 
appropriately, for determining PDF errors later. All of this information is stored in a file 
updated after each Vegas iteration, in case a calculation needs to be stopped or restarted. 

A large number of histograms have been included in the default version of FEWZ. Some 
differential distributions provided include: 

• dilepton transverse momentum; 

• dilepton rapidity; 

• dilepton invariant mass; 

• lepton transverse momenta; 

• lepton pseudorapidities; 

• jet transverse momenta; 

• jet pseudorapidites; 

• AR separation between observable particles; 

• Ht (scalar sum of all transverse momenta); 

• cos(#*), the lepton polar angle in the Collins-Soper frame |14j ; 

• A(f), the lab-frame transverse-plane angular separation between leptons. 

In addition, by reweighting each phase-space point according to certain trigonometric func- 
tions in the Collins-Soper coordinates, we can reconstruct the Collins-Soper moments Aj [T4T 
[T6] , binned in dilepton transverse momentum. We currently support only one-dimensional 
histograms in FEWZ. Since the weights of all events are saved, it is straightforward to extend 
the program to handle higher-dimensional histograms. 



The number of supported PDF sets has been drastically expanded from the previous version 
of FEWZ. All modern distributions are now available. At present, the code supports and 



includes the following sets: 

• ABKM 09 NLO and NNLO [21]; 

• CTEQ versions 6L1 [25], 6.5 [26], 6.6 [27], 10 and 10W [28]; 

• GJR 08 LO/NLO [29J and JR 09 NNLO [30]; 

• MRST 2006 NNLO [31] and MSTW 2008 LO/NLO/NNLO [32 



3.6 PDFs 




All sets but CTEQ 6L1 include PDF error eigenvectors. The results for each eigenvector set 
are calculated and stored in parallel to the central PDF by reweighting each generated phase 
space point appropriately. The results are stored in an auxiliary output file. When the results 
of all sectors are combined, the PDF errors are calculated by summing the eigenvectors in 
quadrature (symmetric sets), finding the standard deviation from the mean (neural network 
sets), or with the procedure described in Ref. [26J (asymmetric sets). 

3.7 Additional Features 

We have provided shell scripts for farming out the sectors in parallel either locally or on 
Condor, and a finishing script which combines the results of individual sectors. In addition 
to the basic operation of combining the sectors and computing PDF errors, the finishing 
script can perform operations such as addition, subtraction, multiplication, and division on 
different runs, all while treating the integration and PDF errors consistently. This can be 
useful for computing the results of disjoint cuts, calculating K-factors, acceptances, and PDF 
correlations between different kinematic regions. These features are described in more detail 
in the manual distributed with FEWZ. 

3.8 Runtime benchmarks 

To demonstrate the performance of the current version of our code, we compare it against 
the previous version by running to various target precisions under typical usage. For the 
following test we use the standard set of cuts described in Sec. HJ minus the isolation cuts that 
were unavailable in the previous version. MSTW 2008 NNLO PDF sets are used, but the 
results are typical for all PDF sets, as PDF evaluation calls require an insignificant amount 
of CPU time. The result is a cross section of approximately 440 pb, leading to an acceptance 
of 46%. These benchmarks are run on an 8-core (2 CPU) Intel Xeon E5335 2.0 GHz machine 
running Scientific Linux 5.5; the new version of the code is able to use the multiple cores 
available on modern hardware. In addition to calculating the cross section with cuts, the new 
version is simultaneously doing all of the histogramming and PDF error eigenset running 
described above; all histograms available have been activated. The integration error versus 
time is shown for both versions below in Fig. [TJ 

For a target precision of 1% and with standard cuts, the previous version of the code 
required 230 hours (9.8 days) and reached a 0.97% relative error. The new version required 
20 hours and reached a 0.84% error. The new version outperforms the old, even on a per-core 
basis, all while computing over a dozen kinematic distributions and the PDF errors for each. 
The detailed results for phenomenological quantities are presented in the next section. 

4 Phenomenological results 

We present in this section phenomenological results for the LHC that illustrate the improve- 
ments in FEWZ detailed in previous sections. Many predictions shown have not previ- 
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FEWZ Error vs. Time 
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Figure 1: Relative error versus time for the previous version of FEWZ and for the current 
version. 

ously been presented in the literature. All numerical results shown use the electroweak 
scheme [34]. The numerical values for the various parameters are shown below: 

M z = 91.1876 GeV, M w = 80.403 GeV, 

T z = 2.4952 GeV, r z _>,+ { - = 0.08399 GeV, G^ = 1.16637 x 10~ 5 GeV" 2 . (7) 

We use 7 TeV for the LHC energy, and use the value of as(M z ) dictated by the appropriate 
PDF set. 

4.1 Benchmark numbers 

We begin by presenting results for the inclusive cross section at NNLO, which we define 
with only the invariant mass cut 66 GeV < Ma < 116 GeV. To define the error arising from 
scale variation, we let the factorization and renormalization scales vary separately in the 
range M z /2 < fj, R>F < 2M Z subject to the restriction 1/2 < jj, R / 'fi F < 2. The PDF error is 
defined using the procedure recommended by the group which produced the fit. We only use 
the most recent PDF sets for which an NNLO extraction is available: MSTW 2008, ABKM 
2009, and JR 2009. We find the following results for the MSTW set, with all sources of error 
indicated: 

MSTW 2008: a inc = 963.7j^|(scale)^;J(PDF) ± 0.5 (tech.) pb. (8) 

We have included as the final error component the technical precision that arises from the 
Vegas integration. This number is below ±0.1%, and does not significantly affect either the 
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central value or the estimate of the error. The PDF uncertainty is the dominant component 
of the error. We now compare this with the central values and errors obtained using the 
other NNLO PDF sets: 

MSTW 2008: a inc = 963.7i|^(PDF) pb; 
ABKM 2009: a inc = 980.5+}^(PDF) pb; 

JR 2009: a inc = 907.3±^(PDF) pb. (9) 

The relative scale and technical errors are similar to the MSTW values, and are negligible 
compared to the PDF uncertainties. We note that the 90% C.L. error bands are shown for the 
MSTW fit, while the ABKM and JR quoted errors should be interpreted as la uncertainties. 
Upon scaling the MSTW uncertainties down by a factor of 1.6, they agree approximately 
with those of the ABKM and JR fits. MSTW and ABKM give similar predictions for the 
central value. The JR 2009 set gives a central value 6% lower than MSTW 2008 and 8% 
lower than ABKM 2009. These results for the inclusive cross section can be compared to 
preliminary results from ATLAS [35] and CMS [36] : 

ATLAS: o inc = 830^(stat.)+gJ(syst.) ± 90(lumi.) pb; 

CMS: a inc = 882^(stat.)+Jg(syst.) ± 97(lumi.) pb. (10) 

The ATLAS result agrees with the theoretical prediction for all three PDF sets within the 
current errors. We note that the CMS cross section is defined with the invariant mass cut 
60 GeV < Mn — 120 GeV, while ATLAS uses the same cut as our default choice. Changing 
the mass window to match the CMS definition leads to a 1.5% increase in our result for all 
three PDF sets; for example, our MSTW cross section becomes 

MSTW 2008: a inc = 977.71^5 ( PDF ) P b > ( n ) 

which agrees with the CMS measurement within errors. 

We now define the standard acceptance cuts, which include the following restrictions in 
addition to the previous invariant mass cut: 

p T ,ie P > 25 GeV, \i]i ep \ < 2.5, 
AR lep>lep > 0.5, AR lep>jet > 0.5. (12) 

We have used the standard definition Ai? 12 = \J (771 — 772) 2 + (4>i ~ 4>2) 2 in defining the 
isolation cuts on the leptons and jets. The partons have been clustered according to the 
anti-A^ algorithm with separation parameter R = 0.5. We note that the simple ARi ep j et 
cut imposed here is not how the isolation requirement is implemented experimentally. This 
effect would also typically be assigned to the efficiency rather than the acceptance. We 
include the cut here to demonstrate the ability of FEWZ to reconstruct jets and constrain 
hadronic activity. The effect of the isolation requirement is to reduce the cross section by 
only a couple of percent, and this cut can be easily removed if desired. The cross sections 



10 



after these cuts have been imposed follows: 

MSTW 2008: a standard = 436.0l}^(PDF) pb; 
ABKM 2009: a standard = 445.6±™(PDF) pb; 

JR 2009: (J s tandard = 404.3^ 1 9 (PDF) pb. (13) 

The relative scale and technical errors are similar to the inclusive case and are much smaller 
than the PDF errors, and have not been included. By taking the ratio of this over the 
inclusive rate, the acceptance can be derived. The scripts provided with FEWZ allow this to 
be done for each error eigenvector to obtain the PDF error. We find the following for each 
set: 

MSTW 2008: A standard = 0.4525i°;°o!o (PDF); 
ABKM 2009: A standard = 0.4544±°;°ol8( PDF ); 

JR 2009: A standard = 0.4456±°X(PDF). (14) 

The technical precisions for each acceptance are well below ±0.1%. The PDF errors on the 
acceptance are much smaller than for the individual cross sections, as expected, and range 
from 0.4% for ABKM 2009 to approximately 0.8% for MSTW 2008. The results for ABKM 
and MSTW are in good agreement, but the JR 2009 acceptance is 1.5% lower than that for 
MSTW, larger than the estimated 90% CL. MSTW PDF error by a factor of two. Scale 
errors for such cuts have previously been studied in Ref. j9], and are at the percent-level or 
below for both the cross section and acceptance. 

To address the limitations imposed by significant phase-space restrictions, we define a 
severe acceptance cut following the analysis in Ref. [T2] : 

p T ,iep > 25 GeV, 1.5 < \rji ep \ < 2.3, 

ARlepJep > °- 5 > ARi epJet > 0.5. (15) 

Only a small slice of the forward region in lepton pseudorapidity is taken. In the study of 
Ref. [1 2j using the old version of FEWZ, an integration precision of only ±3% was obtainable 
after asymptotic running, preventing an accurate estimate of the higher-order corrections in 
this phase-space region. Using the new version of FEWZ, we obtain a technical precision at 
the ±0.5% level after several days of running: 

MSTW 2008: a severe = 37. 09^ (scale) +};^ (PDF) ± 0.18(tech.) pb. (16) 

We note that the difference between the central value presented here and in Ref. [12] is due 
primarily to the 7 TeV energy we use. The technical precision is a factor of a few less than 
the (small) scale dependence, even though the acceptance is only 4% for this cut. We view 
this as evidence that the limitations imposed by integration errors are solved for all studies 
of interest at the LHC The PDF errors are again the dominant uncertainty on this result, 
leading us to study the results for the other NNLO PDF sets: 

MSTW 2008: a severe = 37.09^ (PDF); 
ABKM 2009: a severe = 36.85t£|(PDF); 

JR 2009: a severe = 35.20lg;^PDF. (17) 
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These lead to the following results for the acceptances and associated PDF errors: 



MSTW 2008 



ABKM 2009 



JR 2009 



A 



A 



A 



severe 



severe 



severe 




(18) 



The conclusion of these studies is that FEWZ is capable of providing results with sub-0.5% 
integration errors even in the presence of severe phase-space restrictions, allowing central 
values, scale and PDF errors to be accurately computed for observables of interest at the 



4.2 Results for distributions: inclusive cuts 

We now present results for LHC distributions, both to demonstrate the histogramming fea- 
tures of the code and also to present several new phenomenological results. We begin by 
running FEWZ once for each NNLO PDF set in the inclusive mode, with only a cut on the 
invariant mass 66 GeV < Mu < 116 GeV imposed. Bin-integrated cross sections for several 
standard kinematic distributions of the lepton pair are shown in Fig. [2j Integration errors 
reach a maximum of 1% for each bin, and are much below this value except near phase-space 
edges. Scale errors are small, and the dominant uncertainties come from PDFs. Only these 
are shown in the figure, and are indicated by the hatched bands. The smaller integrated 
result for the JR 2009 distribution is apparent from both plots. The rapidity distribution 
of the reconstructed Z is also slightly flatter for the JR 2009 set. The ratio between the 
first and second bins of the p? distribution differs for each set by an amount larger than the 
estimated uncertainty. However, this bin only includes the range GeV < p?,z <■ 5 GeV, 
and fixed-order perturbation theory is not expected to accurately describe this region. If this 
ratio difference persists after the including the resummation of 1ow-pt,z logarithms, it could 
be an interesting discriminator between different PDF extractions. In addition to distribu- 
tions of dilepton variables, FEWZ also produces histograms of leptonic variables. Results 
for several basic leptonic distributions are shown in Fig. El The Jacobian peak present at 
Pt,i — 45 GeV is visible in the transverse momentum spectrum, as is the tail generated by 
additional QCD radiation beginning at 0{a s ). Since the lepton-pair transverse momentum 
distribution and the lepton transverse momentum distribution above the Jacobian peak are 
generated by additional QCD radiation, their perturbative expansion begins at 0(a s ) and 
FEWZ effectively produces only NLO distributions for these observables. The corrections 
for these quantities could also be obtained from an NLO calculation for Z+l jet, such as 
implemented in MCFM [37]. 

Angular distributions in the Collins-Soper frame [H] yield information on both the cou- 
plings of the Z-boson to leptons, and on the perturbative QCD which produces the Z trans- 
verse momentum. The differential cross section is expressed using the polar and azimuthal 



LHC. 



12 



ft 



100 



80 — 



60 — 



40 — 



20 



-i 1 1 1 1 1 1 1 r 

Inclusive cuts: 

66 GeV < M u < 116 GeV 



"1 1 1 1 1 1 1 1 r 



mi 



mmmmmmmsm 



£2£ 



FTZi MSTW08NNLO 

JR09NNL0 
WZfr ABKM09NNL0 



_i i i i i i i_ 





Y 7 



43 
ft 



300 



250 



200 



150 



100 



50 — 



i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



MSTW08NNLO 
j . . ^ JR09NNLO 

Inclusive M ABKM09NNLO 

66 GeV < M u < 116 GeV 




100 



Pt.z 



Figure 2: Bin-integrated cross sections for the lepton-pair rapidity (upper panel) and trans- 
verse momentum (lower panel) for all three NNLO PDF sets. Only a cut on the invariant 
mass 66 GeV < Mu < 116 GeV has been implemented. The bands indicate the PDF uncer- 
tainties for each set. 
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invariant mass 66 GeV < Mu < 116 GeV has been implemented. The bands indicate the 
PDF uncertainties for each set. 
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decay angles of the lepton in this frame as [T5|[T6 



(1 ° 1 + cos 2 9 + -A (1-3 cos 2 6») + A 1 sin 29 cos . 



dM(i dp^ z dY z dcos 9 d(p 2 

+ sin 2 6> cos 20 + A 3 sin 6> cos + A 4 cos6» + (19) 

where the ellipses denote additional terms ^5,6,7 which we do not consider here. We denote 
the Ai as Collins-Soper (CS) moments. The moments A3 and A4 are proportional to the 
parity-violating couplings of the quark and leptons. A 4 is the only non- vanishing moment at 
LO. All others are generated by additional radiation recoiling against the lepton pair. These 
moments can be obtained using orthogonality relations for the trignometric functions they 
multiply. Defining the moment of a quantity m as 

_ / dcos9d(p m da (Mi h p TiZ ,Y z , cos 9, (p) 
J dcos 9 d(f) da (Ma, Pt,z, Y z , , cos 9, 0) 

we can obtain the CS moments in the following fashion: 

<i ( i-3™«)> = !(*-§). 

1 X 

(sin26>cos0) = -A x , (sin 2 ^ cos 20) = — A 2 , 
5 10 

(sin6»cos0) = -A 3 , (cos0) = ^A 4 . (21) 

We begin by showing the results for the normalized cos 9 distribution in Fig. HI The PDF er- 
rors are completely negligible for this distribution, and all three sets are in perfect agreement. 
This prediction of perturbative QCD is stable against theoretical uncertainties. However, 
we will show later the important effect that acceptance cuts have on the cos# distribution. 
We next display the results for the CS moments in Figs. O and [7J These again show very 
little variation under choice of PDF set or eigenvector within a set. The CS moments are 
also affected dramatically by acceptance cuts, as we demonstrate later. 

In order to reduce backgrounds affecting the Drell-Yan measurement, LHC experimen- 
talists have discussed imposing a cut demanding that the two leptons not be back-to-back 
in the transverse plane [38J. Z-boson events where the leptons are separated by nearly 180° 
suffer from large backgrounds from semi-leptonic b decays. A cut demanding a minimum 
value for this angle reduces such backgrounds. To study the QCD predictions in the presence 
of this cut, we denote by A(f>u the lower cut on the separation angle on the two leptons. For 
example, a A0;/ cut of 3° denotes that all events where the deviation between the leptons 
in the transverse plane is more than 3° are accepted. Two possible problems make this cut 
worrisome from the perspective of stability under QCD corrections. The region A0^ ~ is 
dominated by the emission of soft and collinear gluons, and large logarithms of A0^ inval- 
idate fixed-order predictions for this quantity. Resummation is required. For A0^ > 0, the 
perturbative expansion starts at 0(a s ). Our result is effectively only next-to-leading order 
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Figure 4: Bin-integrated, normalized cos 6 distribution for all three NNLO PDF sets. The 
polar angle of the lepton is defined in the Collins-Soper frame, as indicated by the subscript in 
the plot. Only a cut on the invariant mass 66 GeV < My < 116 GeV has been implemented. 
The bands indicate the PDF uncertainties for each set. We note that the smallness of the 
PDF errors makes it difficult to distinguish the three separate bands in the plot. 



for this observable, and could be obtained by studying Z+l jet at NLO. This indicates that 
the scale variation will be greater than for the inclusive result. 

To study these issues we show in Fig. [8] the results obtained by including all events with 
angular separation greater than some lower cut A<pu at 0(a s ) and 0(a 2 s ); for the inclusive 
cross section these would respectively be the NLO and NNLO cross sections. The onset 
of large logarithms as A<fiu — > is clearly visible in the divergence of the cross sections 
toward the left side of the plot. In this region a resummation of the associated logarithms 
is required. Fixed-order perturbation theory cannot be trusted below a cut value of roughly 
A(pii m 3°, where the bands cross. To study the theoretical uncertainty on the cross section 
and acceptance in the region where fixed-order results should give a reasonable estimate of 
the cross section, we show below the results including scale and PDF errors for the choice 
A(f) U = 3° and with MSTW PDFs: 

^=3° = 903.8lg;f(scale)+|;°(PDF) (pb); 

A A4>=3 » = 0.943^;^(scale)^;S(PDF). (22) 

While the relative PDF errors on the cross section and acceptance are the same as for the 
other cuts studied above, the scale errors on both reach ±2.5%. This is significantly larger 
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Figure 5: Bin-integrated results for the CS moment Aq, presented as a function of the lepton- 
pair transverse momentum, for all three NNLO PDF sets. Only a cut on the invariant mass 
66 GeV < Mu < 116 GeV has been implemented. The bands indicate the PDF uncertainties 
for each set. We note that the smallness of the PDF errors makes it difficult to distinguish 
the three separate bands in the plot. 



than the scale dependence found for the inclusive result or for the standard acceptance 
cut defined previously, and occurs because this obervable begins at one order higher in 
perturbation theory. Although the experimental errors might be reduced by imposing this 
cut, the increased theoretical uncertainty should be accounted for in analyses. 

As an additional phenomenological result relevant to LHC analyses, we study the in- 
variant mass distribution far above the Z-pole. Events in this phase-space region serve as 
a background to searches for high-mass resonances or contact interactions. We run FEWZ 
once, and use the histogramming feature to study the high-mass portion of the Mu distri- 
bution. The result for the range 500 GeV < Mu < 1.5 TeV is shown in Fig. |9l The JR 
2009 PDF set gives consistently lower cross sections than both the MSTW and ABKM sets. 
The MSTW and ABKM fits agree over the entire invariant mass range. In this phase-space 
region, important additional corrections come from electroweak Sudakov logarithms. These 
will be included in a future update of FEWZ. 

4.3 Results for distributions: standard cuts 

To demonstrate the effects of leptonic cuts on the studied distributions, we run FEWZ again 
for each of the three NNLO PDF sets and impose the standard acceptance cuts introduced 
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Figure 6: Bin-integrated results for the CS moments A\ (upper panel) and A 2 (lower panel), 
presented as a function of the lepton-pair transverse momentum, for all three NNLO PDF 
sets. Only a cut on the invariant mass 66 GeV < Mu < 116 GeV has been implemented. 
The bands indicate the PDF uncertainties for each set. We note that the smallness of the 
PDF errors makes it difficult to distinguish the three separate bands in the plot. 
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Figure 7: Bin-integrated results for the CS moments A% (upper panel) and A± (lower panel), 
presented as a function of the lepton-pair transverse momentum, for all three NNLO PDF 
sets. Only a cut on the invariant mass 66 GeV < Mu < 116 GeV has been implemented. 
The bands indicate the PDF uncertainties for each set. We note that the smallness of the 
PDF errors makes it difficult to distinguish the three separate bands in the plot. 
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Figure 8: Results obtained by including a lower cut on the transverse angular separation 
between the two leptons, A</>a, at 0(a s ) and 0(ot 2 s ). Only a cut on the invariant mass 
66 GeV < Mu < 116 GeV has been implemented. 



previously: 



p TMp > 25 GeV, \r]i ep \ < 2.5, 

A.Rl e pj e p > 0.5, ARi e pj e t > 0.5. 



(23) 



We show only a few representative distributions to avoid too severe a proliferation of plots. 
For basic kinematic distributions such as the transverse momenta or rapidities of the leptons 
or Z-boson, the standard acceptance cuts do not dramatically affect the shapes. This is 
demonstrated for the lepton and Z-boson rapidity distributions in Fig. [TUJ The MSTW 
and ABKM Z-boson rapidity distributions become flatter after the standard acceptance cuts 
are implemented, suggesting that their differences from the JR set occur at low-x. No other 
significant differences from the distributions obtained with only the invariant mass cut are 
apparent, besides the obvious changes in distribution endpoints. 

The situation changes drastically when angular distributions of the leptons are studied. 
We begin by showing the cos 9 distribution after standard acceptance cuts are imposed in 
Fig. EJ The shape of the distribution in Fig. H] has completely changed upon addition of 
the cuts. This can be understood by considering the leading-order kinematics. At LO, the 
expressions for the lepton and anti-lepton pseudorapidities are given by 



Vi 



In 



xi 1 + cos 9 
X2 1 — cos 9 



m 



In 



xi 1 — cos 9 
X2 1 + cos 9 



(24) 
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Figure 9: Bin-integrated cross sections as a function of invariant mass for the range 
500 GeV < Mu < 1.5 TeV. No additional cuts have been imposed. 



Events near cos (9 ~ ±1 correspond to events with high values of \rjii\, which are removed by 
the cuts in Eq. ffl2]) . 

The modification of the polar angle distribution, and additionally the lower cut on the 
leptonic transverse momenta, have a dramatic effect on the determination of the CS moments 
defined previously. The orthogonality conditions of Eq. ( 12T]) no longer apply when only a 
finite region of cos 9 is integrated over experimentally. To demonstrate the impact of these 
acceptance cuts on the extraction of the CS moments, we define 'naive' CS moments A™ awe 
which we obtain by integrating over only the region allowed by the standard acceptance cuts 
in the moment definition of Eq. ( 120]) . The result for A™ tve is shown in Fig. [12j Instead of 
beginning at zero for p?,z — GeV and monotonically approaching A ~ 0.75 at px,z — 100 
GeV, A™ lve begins at A™ me m 2, rises to a maximum p T ,z ~ 45 GeV, and falls to A™ me « 1 
at pt,z = 100 GeV. These general features can be confirmed by analytic integration of the LO 
result for the cross section. The acceptance cuts completely change the qualitative features 
of this distribution. A similar effect is obtained for the moment A 2 , shown below in Fig. [13j 
Instead of a distribution monotonically increasing toward 0.75 at large pt,z as seen in Fig. |6l 
it instead monotonically descreases from zero toward A^ aive m —3 in the presence of standard 
acceptance cuts. 
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transverse momentum (lower panel) for all three NNLO PDF sets. The standard acceptance 
cuts of Eq. ffl2|) have been implemented. The bands indicate the PDF uncertainties for each 
set. 
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Figure 11: Bin-integrated, normalized cos# distribution for all three NNLO PDF sets. The 
polar angle of the lepton is defined in the Collins-Soper frame, as indicated by the subscript 
in the plot. The standard acceptance cuts of Eq. ( 1T2|) have been implemented. The bands 
indicate the PDF uncertainties for each set. We note that the smallness of the PDF errors 
makes it difficult to distinguish the three separate bands in the plot. 



5 Conclusions 

In this manuscript we have presented an improved version of the analysis code FEWZ for 
the study of differential lepton-pair production through next-to-next-to leading order in 
perturbative QCD. FEWZ allows the kinematics of the leptons and all associated hadronic 
radiation to be studied, and permits analysis of events both on and off the Z-peak. The 
new features of the code include an efficient, parallelized integration routine that can take 
advantage of distributed computing systems such as Condor. Sub-1% technical precisions 
are easily obtainable even in the presence of severe phase-space cuts that accept only a small 
fraction of the inclusive cross section. Histograms of most interesting kinematic variables 
are now filled during a single run of FEWZ, and PDF errors are automatically computed for 
each histogram bin. Both inclusive results and distributions are obtained with a single run 
of FEWZ on a single multi-core desktop after no more than several days even with extreme 
cuts imposed, or in less time using a Condor system. 

We have presented numerous phenomenological results relevant for LHC studies that also 
demonstrate the new features of FEWZ. We have shown inclusive cross section central values, 
scale variations and PDF errors for all three NNLO PDF fits, and have compared them 
to recent measurements by the ATLAS and CMS collaborations. Additional predictions 
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Figure 12: Bin- integrated results for the naive CS moment iQ me , presented as a function of 
the lepton-pair transverse momentum, for all three NNLO PDF sets. The standard accep- 
tance cuts of Eq. (fl2l) have been implemented. The bands indicate the PDF uncertainties 
for each set. We note that the smallness of the PDF errors makes it difficult to distinguish 
the three separate bands in the plot. 



for representative acceptance cuts are also shown. FEWZ delivers a wealth of additional 
information in the form of differential distributions. We have presented and discussed several 
basic lepton kinematic distributions. 

The most interesting phenomenological results we have obtained are associated with the 
angular distributions of the leptons. We have studied a proposed cut on the transverse- 
plane angular separation between the leptons, which is designed to reduce backgrounds 
in precision studies of the Z-peak. This distribution suffers from large logarithms if the 
leptons are nearly back-to-back, rendering fixed-order perturbation theory unstable. In the 
region where fixed-order calculations can be trusted, imposing this cut increases the scale 
uncertainty from sub-1% to ±2.5%, since the prediction begins at 0(a s ) for this observable. 
The increased theoretical uncertainty must be accounted for in precision studies. We have 
additionally studied the angular distributions of the leptons in the Collins-Soper frame, 
and the associated Collins-Soper moments. These quantities are absolutely stable under 
perturbative corrections, but the effect of acceptance cuts are dramatic, and change their 
qualitative features. 

FEWZ is intended for use in studies of all aspects of lepton-pair production at hadron 
colliders where fixed-order perturbation theory is applicable: for computing inclusive cross 
sections, distributions of basic kinematics and of angular quantities, and especially when 
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Figure 13: Bin- integrated results for the naive CS moment A^ atve , presented as a function of 
the lepton-pair transverse momentum, for all three NNLO PDF sets. The standard accep- 
tance cuts of Eq. (I12p have been implemented. The bands indicate the PDF uncertainties 
for each set. We note that the smallness of the PDF errors makes it difficult to distinguish 
the three separate bands in the plot. 



determined the effects of limited detector acceptance. In order to compare experimental 
results with the most precise available theory, NNLO QCD should be utilized in all stages 
of the analysis, including when acceptance corrections are determined. FEWZ makes this 
possible. It is a flexible framework that makes future inclusion of electroweak corrections and 
other improvements simple. We look forward to its continued use in experimental studies. 
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